Interaction-Induced Magnetization of the Two-Dimensional Electron Gas 
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We consider the contribution of electron-electron interactions to the orbital magnetization of 
a two-dimensional electron gas, focusing on the ballistic limit in the regime of negligible Landau- 
level spacing. This regime can be described by combining diagrammatic perturbation theory with 
semiclassical techniques. At sufficiently low temperatures, the interaction-induced magnetization 
overwhelms the Landau and Pauli contributions. Curiously, the interaction-induced magnetization 
is third-order in the (renormalized) Coulomb interaction. We give a simple interpretation of this 
effect in terms of classical paths using a renormalization argument: a polygon must have at least three 
sides in order to enclose area. To leading order in the renormahzed interaction, the renormalization 
argument gives exactly the same result as the full treatment. 

PACS: 73.20.Dx,71.10.Ca,03.65.Sq,05.30.Fk 



I. INTRODUCTION 



Within the independent-electron picture, the magnetic 
response of a bulk two-dimensional electron gas has two 
sources: Pauli paramagnetism originating from the elec- 
tron spin and Landau diamagnetism originating from the 
orbital electronic motion. After studies of the contri- 
bution of superconducting fluctuations to the magnetic 
response of superconductors above Tc 0:^, Aslamazov 
and Larkin |^ pointed out that electron-electron inter- 
actions make an analogous contribution to the magnetic 
response of normal-metal systems. While the fluctua- 
tion contribution is diamagnetic in superconductors, the 
Coulomb interaction gives a paramagnetic contribution 
to the susceptibility of normal metals; this difference is a 
direct consequence of the different signs of the effective 
interaction in the two cases. 

In their seminal work, Aslamazov and Larkin |^ com- 
puted the interaction contribution to the susceptibility 
of three-dimensional metals and of layered systems at 
zero magnetic fleld. They found that the effect was 
particularly strong for layered systems. In view of the 
importance of the physics of the two-dimensional elec- 
tron gas, the purpose of the present paper is to com- 
pute the interaction-induced magnetization of a strictly 
two-dimensional bulk system. We shall moreover go be- 
yond the zero-field limit considered by Aslamazov and 
Larkin and compute the magnetization for arbitrary clas- 
sically weak magnetic fields. We find that the interaction- 
induced magnetization generally dominates over the Lan- 
dau and Pauli contributions at sufficiently low tempera- 
tures. 

The relevant length scales of the problem are the ther- 
mal length Lt = F / {2ttT) , the magnetic length Lh = 
{h/eBy^^ , the cyclotron radius Rc — meVp/eB, and the 
elastic mean free path iei- Throughout this paper, we 



focus on the regime Lt ^ (ei , which allows us to neglect 
the effects of impurity scattering. Moreover, we restrict 
ourselves to classically weak magnetic fields, defined by 
the condition Lt ^ Rc (or equivalently huJc <C T, where 
LUc is the cyclotron frequency). Within a semiclassical 
approach, this implies that we can neglect the classical 
effects of the magnetic field on the trajectories and only 
need to consider the Aharonov-Bohm phases induced by 
the B field. For classically weak magnetic fields, we can 
distinguish between two magnetic field regimes: The low- 
field regime, considered by Aslamazov and Larkin, where 
Lt ^ Lh and the high-field regime Lt ^ Lh- We 
present analytical results in both regimes and also show 
numerical results bridging these two regions. Most im- 
portantly, we present a simple renormalization argument 
based on classical paths which exactly reproduces the re- 
sult of the complicated full treatment. 

The magnetization of the two-dimensional electron gas 
has been studied experimentally in mesoscopic samples 
|^,|| and in the quantum-Hall regime 1^-^. To the best 
of our knowledge, no experiments have been performed 
on bulk 2d samples at classically weak magnetic fields. 
Such experiments would be a valuable test of our theo- 
retical understanding of the interaction contribution to 
the magnetization. 

This paper is organized as follows. In Sec. || we 
develop the semiclassical approach to the interaction- 
induced magnetization for 2d bulk systems. In Sec. ^ we 
employ the general results derived in Sec. ^ to derive ex- 
plicit expressions for the magnetization in the high- and 
low-field regimes. A curious feature of these results is 
that the interaction-induced susceptibility is third order 
in the (renormalized) interaction strength. In Sec. |rv| , 
we show how the renormalization group approach intro- 
duced in allows one to give a simple semiclassical 
interpretation of this result. We conclude in Sec. ^ by 
comparing the interaction-induced susceptibility to the 



1 



Landau and Pauli susceptibilities and discussing finite- 
size effects. 



II. THE SEMICLASSICAL APPROACH 
A. Basics 

Cooper-channel 

Calculating the interaction contribution to the mag- 
netic response requires one to extend the high-density 
expansion (RPA) of the thermodynamic potential 
by including interaction corrections from diagrams in the 
Cooper-channel. This was first done in the context of 
superconducting fluctuations and then applied to disor- 
dered normal metals P,p^-p^. Such expansions usually 
give reliable results even beyond the high density limit, 
if the relevant sets of terms are properly resummed. The 



relevant Cooper-like diagrams are shown in Fig. 1. The 
screened Coulomb interaction (wavy lines) can be treated 
as local @,||: U{r - r') = AoiV(0)-i(5(r - r'). Here, 
A''(0) = rrif,/ {-kTi^) denotes the full density of states and 
the bookkeeping index Aq = 1 identifies the order of per- 
turbation. For a local interaction, the direct and ex- 
change term are the same up to a factor of (—2) com- 
ing from the spin sums and the different number of 
fermion loops. The straight lines in Fig. 1 represent 
finite-temperature Green functions of the non-interacting 
system. These take the form 

Qr.r'{e) - e{e)G^y{E, + ie) + e{--e)G^^^,{E, + ie) (1) 

in terms of the retarded and advanced Green functions 

The perturbation expansion for this interaction contri- 
bution f2 to the thermodynamic potential, which yields 
the magnetic response, can be formally expressed as 
00 



^ — ^ dri. . .dr„ i;^(ri,r2) . . . Ej^(r„,ri) 

i^Tr{ln[l + AoE.]} . 



(2) 
(3) 



Here, u denotes the bosonic Matsubara frequencies w = 
2TTm/P (m is any positive or negative integer) with 
P = l/ksT. (We employ units such that fc^ = 1 in this 
paper.) The particle-particle propagator Tt^ is expressed 
(in position representation) in terms of products of finite- 
temperature Green functions as pT[ 



(4) 



where the sum runs over the fcrmionic Matsubara fre- 
quencies e = (2n -I- l)7r//3. The short-length (high- 
frequency) behavior is included in the screened interac- 
tion, thus requiring a cutoff of the frequency sums at the 
Fermi energy uM. 



Semiclassical Green function 

In view of the fact that the Fermi wavelength is the 
smallest length scale in the problem, our strategy will be 
to replace the free Green functions in the particle-particle 
propagator by their semiclassical approximations. Gen- 
erally, the semiclassical approximation to the Green func- 
tion Gr,r'{Ep+ie, B) is expressed as a sum over all clas- 
sical paths from r to r' at energy Ep mm. For the bulk 



2d electron gas, there is only a single such trajectory, 
namely the straight line connecting the two points. For 
e = and B = Q one therefore has 



1 me e-xjp[iSr^r'{EF)/h — «7r/4] 



/2^ h-" 



(fcF|r-r'|)i/2 



(5) 



where Sr.r'iEp) = hkp\r — r'\ is the classical action along 
the path. Moreover, since we assume the magnetic field 
is classically weak, the field affects the action along the 
path through 



SiB) = S{B = 0) + -^j^ dx A(x) 



(6) 



where the integral is along the unperturbed straight line 
path. Finally, it turns out that only small values of the 
imaginary part of the energy e should be considered, so 
we can use the relation 



{dS/dE) = t 



(7) 



where t = \y — y'\/vp is the time of flight from r to r'. 
In this way, one obtains the semiclassical Green function 
for finite field and finite (imaginary) energy e: 



G^,^,{Ep+ie,B) = G^^,{Ep,B^Q)eyiY> 



he 



dx A(x) 



exp 



hvf 



(8) 
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Semiclassical particle-particle propagator 

In the calculation of the thermodynamic potential Eq. 
one should neglect all rapidly oscillating contribu- 
tions in I];j(r, r') as these will give a small contribution 
upon integration. Thus, in the particle-particle propaga- 
tor S, it is necessary to pair advanced and retarded Green 
functions, and furthermore to pair each path in the semi- 
classical expression for with those in G"^ for which 
the dynamical phase factor cancels. The obvious case of 
pairing each path with itself is excluded because it yields 
no field dependence in Q and hence zero magnetization. 
Thus one is led to consider pairs of time-reversed paths- 
for these the dynamical phase cancels but the magnetic 
field part is multiplied by two. The pairing of with 
G^ means concretely that one should keep only those 
terms in which e and lu — e have opposite sign in the sum 
over Matsubara frequencies. Using the relation 



e{uj-e)<0 



\2e~Lj\t/h] 



exp[-~\uj\t/h] 
smh(2tTT / pn) 



(9) 



f(R) 



exp I -R • (p 



2e 



A) 



(14) 



and therefore with its generators 

- {Px - (2e/c)A^) 
= iPy - i'^e/c)Ay) . 

Noting that, first, S„(r,r') is invariant under rotation 
and, second, the Landau Hamiltonian for a particle of 
charge (— 2e) can be written as 



Hl 



Ul+fll + ^l, (15) 



one obtains the final result for S 



we see that Ea;(r,r') is diagonal in the basis {V'n,;} of 
the eigenvectors of Hl and Jz, where n and I are the 
Landau level and angular momentum quantum numbers, 
respectively. For / = the Landau level wavefunction (for 
a particle of charge — 2e) has the well-known form 



E^(r,r') = SI?'(|r- 



, J 2ie 



dxA(x)l (10) 



where we have introduced the zero field limit of the 
particle-particle propagator 
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r - r' = 



exp[— |m||r — t'I/Lt] 



47rLT|r-r'| sinh(|r - t'I/Lt) 



B. Derivation of the eigenvalues 



The ladder-diagram contribution to the thermody- 
namic potential is expressed in Eq. (^) solely in terms 
of the eigenvalues of the operator 'S^{r, r'). We therefore 
need to solve the eigenvalue equation 



dr' I]^(r,r')Vn,i(r') = CVn,i(r) 



(12) 



where n and / are quantum numbers. 

Assuming from now on the symmetric gauge A = 
B X r/2, Eq. ^) reads 

S4r,r')-sL"H|r-r'|)exp|^|(rxr')} . (13) 

It can be easily checked that any operator of the form 
Eq. (^|) commutes with any element of the magnetic 
translation group 



^n,o(r) = exp 



2Ll 



r2 



(16) 



with Ln the Laguerre polynomial and Lh — {fi/ eBY^'^ . 

Finally, an important property of the r(R) is that they 
commute with both Hl and Stj(r,r') but not with J^. 
Since, moreover, within a Landau level there is no stable 
(11) subspace for all the r(R), the eigenvalues cr"'' cannot 
depend on the angular momentum quantum number I. 
At r = 0, Eq. (pi), then, reads 



1 



drSW(|r|)^„,o(r 



(17) 



Using the explicit expressions for y!"^\v) [Eq. (^t|)] and 
^nfl, we finally obtain for the eigenvalues 



n,l 



(18) 



" dx ^""P^ '""'^^ exp{-a:V2a^}L„(xVa') 
sinh X 



where a = Lh/Lt is the essential dimensionless param- 
eter. It is important to keep in mind that the screened 
interaction already implicitly takes into account the ef- 
fect of the interaction on scales shorter than the Fermi 
wavelength, so the integral over x should be cut off for 
small X at approximately Xmin = '^/{kpLT)- 
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C. Reordering of the sum 

The interaction contribution to the thermodynamic 
potential is given in terms of the eigenvalues cr" by 



(19) 



Here we have already taken proper account of the degen- 
eracy of the eigenvalues by the prefactor 2BA/(f>Q. The 
magnetization per unit area now follows by differentia- 
tion with respect to B, 



ln(l + Ao< 



1 + Xo<JJ\ dB 



(20) 



When done naively, the sum over the quantum number 
in this expression diverges. We assume that this is asso- 
ciated with the inadequate treatment of the interaction 
at short distances. We expect that when working with 
the full screened interaction, the contribution of large 
quantum numbers is appropriately suppressed. Hence, 
we reorder the sum in such a way that the sum becomes 
convergent and the eigenvalues with sufficiently large n 
do not contribute appreciably to the sum. This philoso- 
phy is completely analogous to the approach taken in the 
work on the fluctuation contribution to the diamagnetic 
susceptibility in superconductors above Tc In fact, 

our reordering closely follows the reordering proposed by 
Payne and Lee 0]. 

In a first step, we compute 



L'^ieBpyn) - -Lr^ieBp'/n) 



(21) 



To simplify this expression, we use the recursion relations for Laguerre polynomials 



xL'^ix) = nLn{x) - nLn-i{x) 

xL'^{x) = {n+ 1)L„+i(.t) - (n -I- 1 - a;)L„(a;) 



(22) 
(23) 



and obtain 



2i3^ = {n + 1)K ■'-<y-.\+ nya-^ - 



(24) 



We can also rearrange 



^ (n + 1) In /|ti = ^ „ In /„ _ ^ („ + 1) In /„ 

n=0 n=0 n=0 



n=0 



(25) 



Using these expressions, we have for the magnetization 



M ^ 11 V V(n + 1) |2 m ^ + - ^°[-""^ - - } , 



(26) 



In terms of the notation 



X, 



1 + AofT" 



(27) 



we have 



(28) 



For all cases considered below, X" <C 1 so that 



(29) 
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This expression will be our starting point for computing 
the magnetization and the susceptibility. Referring back 
to the definition of X" in Eq. ( p7| ) above, we see that all 
contributions to the susceptibility are at least third-order 
in the interaction Xq. 

III. MAGNETIC SUSCEPTIBILITY 

Using the general results of the last section, we now 
find expressions for the susceptibility in two limits- small 
and large magnetic field- and then evaluate the suscep- 



tibility in the intermediate regime numerically. Before 
considering the various regimes, it is useful to note that 
X", as defined in Eq. (p7|), consists of two factors with 
noticeably different behavior. On the one hand, because 
the integrand in Eq. ( p^ ) behaves as l/x at small x, 
both cr" and Ao/(l + Aqct") are dominated by a loga- 
rithmic singularity at zero and so have little magnetic 
field dependence. On the other hand, using the relation 
Ln{x) — Ln-i{x) = —xL\_i{x) jn {L\ is a generalized 
Laguerre polynomial), we can rewrite Act" = cr"^^ — cr^ 
as 



1 1 ^oo p-\m\x 

AaZ - -7 -TTT^ / dxx''^—e^p{-xy2a^}LlXx''/a^). 



Here the x integration is well behaved at small x, and so For ri'^nQ — c? ■ max(l, |m|), Eqs. ( p2|) and (33) yield 
the lower limit Xmin can be replaced by zero. |l8| 



A. Small-magnetic-field (high-temperature) limit 

The small-magnetic field, or equivalently high- 
temperature, limit is defined by a ^ 1. The factor 
^-\ra\x j gjj^jjj. provides an upper cutoff at min(l, |m|~^) 
in the integrals Eqs. ([l^ ) and (|30|). In addition, x^ ja^ is 
much smaller than one in the entire range of integration. 
We can therefore use the asymptotic expression [0 

e-'^I^KXx) . 



-^(^.x/4)'"/V„(M'/'), (31) 



valid in the range 0<a;<n^/'^ \Ja[x) denotes the Bessel 
functions, r(n) the Gamma function, and v = 4n+2a+2] 
and obtain 



sinh X 



Jo (2^n+ l/2a;/a) (32) 



1 



1 



Vn + 1 2a 



dxx- 



\m\x^ f 2\/n+T X 



sinh: 



-Ji 



I /•niin(l,|m| ) ^ 

- — = - ln(fcFiT/niax{l, |to|}) 

2 . L, . X 2 



Aa" 



1 

2^ 



dx X 



-|m|a 



sinh: 



(34) 



(35) 



up to constants of order one, and so the n dependence can 
be neglected. For n>no, both Acr^ and depend on 
n. However, for cr" the dependence is only logarithmic, 
since it merely amounts to replacing the upper bound of 
the integral by a/y/n. Hence the dominant n dependence 
of comes from Ac". 

From these results for the eigenvalues, the magnetiza- 
tion [Eq. (29)] to lowest order in the small parameter 



(33) In (kpLT) is 



OO 

E 



dx 



xexp(— |m|a;) 
sinh X 



J I {2^/nx/a) 



(36) 



The sum over n converges only slowly and of order no 
terms contribute. In view of the fact that n is multiplied 
by x^ I o? in the argument of the Bessel function, we can 
replace the sum over n by an integral. This yields the 
final expression 



M = 



Trln (kpLx) 



(/cF^t) IXLandaul^ 



(37) 



where |xLandau| = e2/127rmc^ and Ct is given by 

OO /.OO J 



where we define 

f^{n) = ^/n 



' xexp{-\m\x) , I- 

dx — J\ (Zynx 

smh X 



(38) 



(39) 
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We see that the interaction-induced contribution to the 
magnetization will generally be larger than the Landau 
magnetization due to the large factor kpLr The fac- 
tor l/ln(fcFir) must be interpreted as a renormalized 
interaction strength in the Cooper channel P,p^. It is 
interesting that the interaction contribution to the sus- 
ceptibility is third order in this renormalized interaction 
strength. This unusual state of affairs can easily be un- 
derstood by considering the classical paths involved, as 
we shall discuss in Section 

The magnetic susceptibility % obtained above is at zero 
field and to lowest order in the renormalized interaction 
strength which is proportional to 1/ InikpLT)- It is possi- 
ble to derive, with a similar approach, an expression for 
x(0) without expanding in the renomalized interaction 
strength. This is done in Appendix ^ one obtains 



X(0) 



I XLandau 



where we have introduced 



An 



2 + Ao(2a"(r)-.9^(0) 



\rn\x 



2a" = 2ct"(B = 0) = / dx^ , 

Jx^i„ sinha; 

~ hi{kpLT/ max(l, |m|)) . 



(40) 

(41) 
(42) 

(43) 



B. Large-magnetic-field (low-temperature) limit 

In the high-magnetic-field, or equivalently the low- 
temperature, limit defined by a ^ 1, the factor 
exp{— x^/2a^}L„(a;^/a^) always cuts off the integral in 
Eq. (|8|) at X't^l so that we can approximate sinhx ~ x. 
Hence, we find 



I |-min[|m| Sq/vH 

- nim[ln{kpLH /\/n), \n{kpLT/\m\ 



(44) 



2(n+ 1) 



rfj^ye-Hm|v-vV2^1(y2). (45) 



In this case, the sum over n converges rapidly (faster 
than 1/n^), but typically about 1/a terms contribute to 
the Matsubara sum. Neglecting again the logarithmic 
dependence of cr" on n and \m\, we can make progress 
by noting that the sum over m can be turned into an 
integral. This yields the high-field result 



M = 



C 



H 



TT \v?(kpL}i) 

where the constant Ch is 



(fcF-^_ff)|XLandau|S, 



(46) 



/OO ^ 
dm 
-OO n 



^ (n-M)2 



0.74 



(47) 



with 



F„(m)= / dyye-\^\y-y''^L\{y^). 



(48) 



The principal difference between the results for high and 
low fields is thus the replacement of the thermal length 
Lt by the magnetic length L}j. 



C. Intermediate range 



When a is neither much smaller nor much larger than 
one, it is not possible to obtain a simple expression for the 
magnetic response. In this regime, we have performed a 
numerical integration of Eq. (|l^) to compute the eigen- 
values cr" as well as their derivatives with respect to B. 
The magnetic susceptibility is then obtained through the 
field derivative of Eq. 



Rrh^ X! ' 



{XZ-^fdXZ-^/dB 

(i + xZ]-i)2 



(49) 



by direct summation over the eigenvalue index and the 
bosonic Matsubara frequencies. 

Fig. H shows the resulting x/|xLandau| as a function 
of magnetic field h = = {2TrL'^/<j)Q)B for fixed val- 
ues of the temperature. The cross-over between the low 
and high field regimes is clearly seen. As B increases, 
X has a slight maximum around Lh ~ Lb which arises 
from the competition between the increased field sensi- 
tivity of large triangles and the thermal suppression of 
long sides. In the large-field regime, the numerical re- 
sult is in reasonable agreement with the value obtained 
from the asymptotic expression (^6|): for Lt — 10Lh and 
kpLn — 64:, x/|xLandau| ~ 0.28 numerically and 0.18 ana- 
lytically. 

Fig. H shows the temperature dependence of x at fixed 
Lh/Lt in the low-field regime. The (kpLx)/ I'n'^ (kpLx) 
behavior is apparent, particularly in the inset. Again, 
the numerical result agrees nicely with the asymptotic 
result: at the low temperature T/Ep — 10^^, Eq. ( |37| ) 
yields x/IXLandau| ~ 0.51 while our numerical result is 
0.52. 
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IV. SEMICLASSICAL INTERPRETATION 

We have seen above that the magnetic response has a 
rather pecuhar property: it is third- order in the renor- 
mahzed couphng constant A — 2ln''^ {kpLT^n), both in 
the low and high field regimes. Within the approach used 
up to now, it is difficult to understand the physical ori- 
gin of this behavior. In this section, we show that an 
approach in terms of classical paths provides a natural 
understanding of this fact. Considering for instance the 
low field regime (the argument can be transposed to high 
fields with no essential difficulty), we shall moreover re- 



cover precisely the expression Eq. (|37|) in a much simpler 
way. 

From the expression Eq. (|l^) for the particle-particle 
propagator, it is clear that the interaction contribution 
to the thermodynamic potential in Eq. can be writ- 
ten as a sum over closed polygonal paths, where each 
vertex is associated with an interaction event and the 
magnetic field enters only via the Aharonov-Bohm fac- 
tor associated with the magnetic flux enclosed by the 
polygon. Performing the sum over bosonic Matsubara 
frequencies in Eq. (j^) and grouping the field dependent 
terms together, one finds 



(50) 



n=l 



(3 n 



dri . . . dvn E(ri, r2) . . . E(r„, ri) coth 



itot(ri 



2Ln 



47rAtot(ri, . 



^)B 



(51) 



where E(r, r') is the particle-particle propagator for B = and u; = Q, itot is the total length, and At^t the area 
enclosed by the polygon. 

As it stands, Eq. (|5^) is of little practical use because the series in n is strongly divergent (the term of order n 
is typically larger than the one of order n — 1). It is, however, possible to apply a simple renormalization-group 
argument, introduced in [0 and discussed in more detail in H^. Indeed, as we already stressed when deriving the 
eigenvalues of the particle-particle propagator, E(r, r') must be cutoff at Aq — \/kp because the use of the screened 
interaction assumes that all high-momentum degrees of freedom have already been integrated out jl^ . We therefore 
have 



E(r,r';Ao) = — — .\ 777^ for|r-r'|>Ao 

A'kLt]^ — r I smh(|r — r |/Lt) 

= for |r-r'| < Ao . 



(52) 



Refs. show that if a new length scale A > Aq much smaller that any other characteristic length scale of the 

problem [Lt or Lh) is introduced, one can replace Aq by A in Eq. (52) provided the "bare" coupling constant Ao in 
Eq. (pG) is replaced by the renormalized one 



AflG(A) 



Ao 



1 + (Ao/2) ln(A/Ao) ' 



(53) 



For the leading behavior in \Dr^{kpLT), we can assume A = cLt with e ^ 1 but assumed fixed as \iY{kpLT) goes to 
infinity. In that case AflG(A) = 2 \nr'^{kpLT) + 0(ln^^(fcj?LT)) is small, and Eq. (|50| ) becomes a genuine perturbative 
expansion whose leading behavior is given by the first non- vanishing term. Clearly, f2*^^^ = 0, and O^^-* is independent 
of the magnetic field and so does not contribute to the magnetic response. Therefore, the leading behavior is given 



by 17^^^- that is, third-order in In ^{kpLr)- as illustrated in Fig. |[ and we have 



lii(A; iJ-Lr) 



L2 ,L23 >A 



dradrg S(0, r2)E(r2, r3)S(r2, 0) coth 



mi 



(54) 



{L2 = |r2|, L3 = ksl, L23 = \r3 - r2|). Expressing aU 
distances in the integral in units of Lt gives 

Crg 



C 



RG 



xiB = 0) = 



with the constant Crg given by 



XLanda 



(55) 



dr2dr3 



1 



Al^ coth(Ltot/2) 



L2L23L3 sinhL2 sinhi23 sinhLs 



(56) 



Because of the factor A? ^ in the numerator, the inte- 
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grand here is regular and the cutoff can be taken to zero. 
By considering the limit Aq ^ ln(/cyLr) ^ 1 (instead 
of Ao = 1) where both the standard and renormaliza- 
tion group approaches are accurate, it can be shown that 
Cbg = Ct, so that Eq. ( |55| ) is strictly equivalent to 
Eq. ( p7| ) . This approach shows clearly that the third power 
of the coupling constant arises because only trajectories 
with three or more vertices enclose flux. 

V. DISCUSSION AND SUMMARY 

In this paper, we study the interaction contribution to 
the magnetization of a two-dimensional electron gas in 
the limit Lt ^ Rc,ic\- We find that this interaction- 
induced contribution is paramagnetic for the repulsive 
Coulomb interaction and dominates over the Landau 
diamagnetism at small enough temperatures and fields 
. The Pauli paramagnetism is even smaller than Lan- 
dau diamagnetism in GaAs/AlGaAs heterostructures be- 
cause of the small effective mass and the reduction in the 
g factor. It appears from the quantitative answer that 
one needs to go to rather low temperatures before the in- 
teraction becomes truly larger than the Landau suscepti- 
bility (cf. Fig. H). Still, such temperatures are possible in 
two-dimensional electron gas systems. It should be pos- 
sible to distinguish the interaction contribution by way 
of either its temperature dependence (since the Landau 
susceptibility is T independent) or its dependence on kp 
in a gated structure. 

We find that the leading contribution to the magne- 
tization is third-order in the renormalized Coulomb in- 
teraction. This can be given a natural semiclassical in- 
terpretation in terms of the classical-path picture for 
the thermodynamic potential: a polygon must have at 
least three sides in order to enclose area. Moreover, the 
much simpler classical-path approach gives exactly the 
same answer as the eigenvalue calculation. Higher-than- 
third order contributions (in the bare coupling constant) 
predominantly lead to an substantial downward renor- 
malization of the coupling constant. This picture has 
been made precise in the present paper by means of a 
renormalization-group approach. 

In the low- field (high-temperature) limit Lt ^ Lh, 
the temperature dependence of the susceptibility is 
l/Tln^T: this comes from the thermal length Lt which 
dominates here because only trajectories shorter than Lt 
make significant contributions to the Green functions. 
At the lowest temperatures, this behavior is cut off by 
a finite magnetic field once Lh < Lt- In the high- field 
(low-tcmperature) limit Lt S> Lh, the susceptibility is 
no longer temperature dependent. The Green function is 
still dominated by trajectories shorter than Lt, but now 
trajectories enclosing more area than L^ contribute with 
random signs due to the Aharonov-Bohm phases. Hence, 
in this case, the relevant cutoff length is Lh- 



So far, we ignored dephasing due to inelastic scatter- 
ing. At low temperatures, this should be mostly due to 
electron-electron scattering. We expect however that de- 
phasing will not significantly affect our results. Within 
the semiclassical approach employed in this paper, de- 
phasing suppresses the contribution of trajectories longer 
than the dephasing length L^. For a clean Fermi liquid 
such as discussed here, one expects L^ ~ Thus, at 

sufficiently low temperatures the dephasing length should 
always be longer than the thermal length Lt ~ 1/T. 
Correspondingly, the suppression of trajectories due to 
thermal smearing should always set in before the sup- 
pression due to dephasing. 

It is interesting to compare the present results with the 
contributions to the susceptibility of (chaotic) mesoscopic 
samples of linear size L within the independent-electron 
approximation ||l9|-|2l| 

I, (57) 

and due to interactions 

(kpL) 

X ^ IXLandaul, (58) 

where we have taken Lt^L, T> A (A is the level spac- 
ing), and Lh^L. In contrast to the bulk results derived 
in this paper, these expressions are zeroth or first order 
in the renormalized interaction constant. These contri- 
butions exist for mesoscopic samples because, in finite- 
size systems, flux-enclosing trajectories are produced by 
scattering from the geometric boundaries of the system. 
Nevertheless, apart from the different order in the renor- 
malized interaction, the finite-size result due to interac- 
tions is qualitatively the same as the bulk result derived 
here. 

Orbital magnetism in mesoscopic samples has been a 
controversial issue over the last decade, both for bal- 
listic and diffusive structures, ring and dot geometries 
P,p|j2^[2^. In particular, the fact that the measured 
values are apparently substantially larger than the theo- 
retical results has attracted a lot of attention. In order 
to benchmark the theory in a simpler system, we think 
the magnetization of a clean two-dimensional electron gas 
should be measured and that this would provide valuable 
information in addressing the "persistent current prob- 
lem" in rings and dots. 
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APPENDIX A: 

In this appendix, we derive an expression for the sus- 
ceptibihty at zero field, Eq. (|4^), which does not involve 
expanding in the renormalized coupling constant. Start- 
ing from the magnetization Eq. (p9[), the essential in- 
gredient needed is expressions for the eigenvalues. One 
can check that both the Bessel approximations Eqs. (^2|)- 
( ^ ) and the replacement of the discrete sum over Landau 
level index n by an integral only yield corrections of or- 
der to the magnetic susceptibility. Therefore, as long 
as we are only interested in the susceptibility at B = 0, 
we can make the change of variables n ^ = bn where 
b = = (27rL^/(/)o)i3 is proportional to the magnetic 
field. We can thus write 



2Act: 



ri-l 



(Al) 



where /^^(C) is defined by Eq. (|39|). In the same way, 
taking the derivative of Eq. ( p2[ ) with respect to b yields 



--i/((n- 1/2)6) ^-i[/(C)+0(6)]. (A2) 



Thus, using the notation of Eqs. (^) and ( p3|) 



2cr,. 



2ct° - g{(n - 1/2)6) ~ 2a" - g{0 + 0{b) 



(A3) 



From these expressions for the eigenvalues, we obtain 

Ao/^(0/e 



x::-' = xu^,b) = -b 



2 + Ao(2a"-.g(0) 
bXM)U0/^ + O{b^ 



(A4) 



Note that it is necessary to compensate the factor b~^ 
originating from the term ndn when changing variables 
from n to ^. Using the above expression, one can there- 
fore write the magnetization up to corrections of order 
b\ 



M 



— E 



to 

which immediately gives Eq. (Eol 



X^iC,bf + 0{b^) 



(A5) 
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FIG. 1. Leading Cooper-channel diagrams for the interaction contribution to the thermodynamic potential. Because we 
can take the interaction to be local (a 5- function), the corresponding Fock-like diagrams differ from the Hartree-like diagrams 
shown only by a factor of —1/2. 
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FIG. 2. Magnetic field dependence of the susceptibility at two temperatures, _Ef/T = 2000 (solid) and 200 (dashed). The 
Landau susceptibility xl is the natural unit for x; the lower panel shows the same data on an logarythmic scale. Note the 
cross-over in behavior when Lt^Lh- 
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FIG. 3. Temperature dependence of the interaction-induced low-field susceptibility; Lh/Lt = 8.9 is fixed. The total 
susceptibility differs from the interaction-induced contribution shown here only by a constant offset due to the Landau and 
Pauli contributions. The inset shows the same data on a log- log scale. For this low field, the contribution of the w = 
Matsubara frequency (dashed) gives a substantial portion of the result (solid). Note the approximate power law increase in x 
at low temperatures, consistent with the asymptotic expression l/Tln^T. 




(a) (h) 

FIG. 4. (a) Typical trajectory of lowest order in the coupling constant which contributes to the interaction contribution 
to the magnetization. Each vertex of the polygon corresponds to an interaction event. Note that at least three interaction 
events are required to obtain a trajectory enclosing magnetic flux. This explains the fact that the interaction contribution 
to the magnetization is third-order in the (renormalized) interaction, (b) Higher-order contributions predominantly lead to a 
renormalization of the third-order result due to short trajectories like the one shown. 
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